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As a qubit usually has a limited lifetime, its manipulation should be as fast as possible, and 
thus non-adiabatic operation is more preferable. Moreover, as a qubit inevitably interacts with its 
surrounding environment, robust operations are of great significance. Here, we propose a scheme for 
quantum manipulation of the polariton qubit in circuit QED using non-adiabatic holonomy, which is 
inherently fast and robust. In particularly, the polariton qubit is shown to be robust against arbitrary 
low-frequency noise due to its near symmetric spectrum, which can also be convenient manipulated 
by external microwave driven fields in a holonomic way. Therefore, our scheme presents a promising 
way of manipulating polariton qubits for on-chip solid-state quantum computation. 


I. INTRODUCTION 

Geometric phases depend only on certain global ge¬ 
ometric property of the evolution path, and thus are 
largely robust against local noises. This distinct fea¬ 
ture makes them promising in implementing quantum 
computation in a fault-tolerant way [lj and high-fidelity 
geometric quantum gates have been experimentally re¬ 
alized UM- Geometric phases can be classified by two 
types, i.e., the abelian and non-abelian ones, quantum 
computation based on which are called geometric JHzl 
and holonomic quantum computation (HQC) [ffl , re¬ 
spectively. The implementation of quantum gates with 
geometric phases can be achieved by both adiabatic and 
non-adiabatic evolution. However, the main challenge 
of the adiabatic method is the long run time, which 
is comparable with the lifetime of typical qubits Si. 
For this reason, non-adiabatic geometric gates should be 
more preferable. Recently, non-adiabatic HQC has been 
proposed with three-level lambda systems |12l4l5| with 
several experimental verifications of elementary gates for 
universal quantum computation [l6l - [l9j . However, in 
such schemes, the upper excited state is resonantly cou¬ 
pled, and thus the limited lifetime of this state is the 
main challenge in practical experiments. 

Here, we propose to manipulate the polariton qubit 
in a non-adiabatic holonomic way. The polariton qubit 
is choose to be a subset of a V-configuration three-level 
system, which is formed by three dressed states in a typ¬ 
ical circuit QED setup. Moreover, due to the polari- 
tonic nature of the qubit, the polariton qubit is shown 
to be robust against arbitrary low-frequency noise. Fi¬ 
nally, the non-adiabatic holonomic quantum gate can be 
constructed by external driven microwave fields. There¬ 
fore, our scheme presents a non-adiabatic and geometric 
way of manipulating polariton qubits for robust on-chip 
quantum computation. 
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II. THE POLARITON QUBIT 

We firstly consider implementing the pro pos ed polari¬ 
ton qubit with a typical circuit QED setup [20(, where a 
superconducting transmon qubit [2l| is capacitively cou¬ 
pled to a one-dimensional transmission line resonator (ID 
cavity). Setting h = 1, the interaction Hamiltonian can 
be written as [22| 

Hjc = w a |l)(l| + w r a) a + g(aa + + (1) 

where w a and w r are the frequencies of qubit and cav¬ 
ity, respectively; a~ = |0)(1| with |0) and |1) being the 
ground and excited states of the transmon qubit. Label 
|n) r as the Fock state of the cavity, the ground state of 
Hamiltonian dH) is |G) = |0)|0) r = |0, 0) with its eigen¬ 
value set to be Eg = 0. As shown in Fig. la, for n > 1, 
due to the transmon qubit-cavity interaction, each energy 
level is spitted into two, the eigenstates are 

| —, n) = cosa„|0, n) — sina n |l, n — 1), (2a) 

|+, n) = sino:„|0, n) + cosa„|l, n — 1), (2b) 

and corresponding eigenvalues are 

E n ,± = mo? + i ± \/5 2 + 4 ng 2 ^ , (3) 

where tan(2a n ) = 2g^/n/5 with the detuning 6 = w r —w a . 
The transition frequencies are 

w n ,± = w r ± 7}(y/5 2 + 4 (n + 1 )g 2 - y/6 2 + 4?ig 2 ),(4a) 

=UJ v + \{\]S 2 + 4 (n + 1 )g 2 + y/6 2 + Mrg 2 )(4b) 

U} n ,\ = Ur - + 4(n + 1)<? 2 + y/8 2 + 4np 2 ),(4c) 

where are the transition frequencies between the 

same transmon states, while co n y and \ corresponds 
to the transitions with different transmon states. It is 
well known that the dynamics of Hamiltonian ([T]) is con¬ 
strained within the subspace with certain |n) r . For the 
zero- and one-excitation subspace, the three lowest eigen¬ 
states |G), | —, 1), and | + , 1) form a V-configuration, and 
the latter two are defined as the polariton qubit, as shown 
in Fig. 1. For the sake of simplicity, we write the qubit 
states |±, 1) as |±) in the following. 
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III. DECOHERENCE 


The major decoherence source in superconducting 
qubits is the low-frequency noise [2314251 ] . The three- 
level system we considered has better stability under 
low-frequency noises compared with conventional super¬ 
conducting qubits. We first consider the influence of 
the transverse noise H' = A x af from a superconducting 
qubit on the polariton qubit. The low-frequency nature 
of the noise determines that it cannot resonantly cou¬ 
plings the transitions |±) O |G) of the polariton qubit. 
Moreover, (±|i7'|q=} = 0, i.e., it also can not induce qubit 
states transition. Therefore, the noise can be treated 
as static fluctuations with perturbation theory and the 
corrected eigenenergies of the qubit states, up to second 
order, are 


El,± = E h± + H[ >± + En Ea e^-eI - (5) 

Note that H[ ± . n a = 0 for n > 2, and thus will greatly 
simplified the caculation of Eq. ©. We rewrite E\ t ± = 
Ei t ± +SE\ t ±, and the correction energy of E\_ and E\,+ 
is calculated to be 


SE h- = 

^El,+ — A x ( 2uj 0i/ / 4wi i+ _ 4<v 1|X )■ 


(6a) 

(6b) 


They can be represented by u> r and g as 

SEi,- = ^ 

5Ei t+ = 


2 

At 


_1_1_ 

uiT-g u) r +g-2g 2 /(u r +g) 


ui T +g ui r -g-2g 2 /(uir—g) 


(7a) 

(7b) 


As j<Sw r , they can be approximated as 
5E X ± » T A 2 x g/ul, 


( 8 ) 


and thus the correction of the qubit energy splitting is 
5E* = 2A x g/ujl. Qubit dephasing is determined by this 
correction and is significantly suppressed by a factor of 
~ A x /uj%, similar to the case of conventional supercon¬ 
ducting qubits working at optimal points [26ll28[ . How¬ 
ever, the optimal points are not immune to longitudinal 
low-frequency noise, which couples to the qubit through 
diagonal matrix elements and thus generates a shift of 
the energy splitting. Fortunately, projecting longitudinal 
low-frequency noise of the superconducting qubit into our 
polariton qubit space, it will also be transformed to the 
transverse type of noise, and thus has minor contribution 
to the decoherence of our polariton qubit, as the case dis¬ 
cussed in the transverse noise. Therefore, the polariton 
qubit considered here is protected from both transverse 
and longitudinal noises. This is in stark contrast to that 
of conventional superconducting qubits where the first- 
order perturbation always plays the leading role. This 
is ensured by polaritonic nature of the qubit as well as 
the symmetric energy spectrum, where although the per¬ 
turbations may induce the coupling between qubit states 
and other energy levels, their contribution to the splitting 
energy are almost cancelled. 



FIG. 1. The polariton qubit states, (a) Level structure of 
the qubit-cavity coupled system, (b) Under approximate pa¬ 
rameters of the driven field, the interaction of the lowest V- 
configuration three levels can be used to realize single-qubit 
gates, (c) The effective interaction in the eigenbasis of the 
V-configuration levels. 


IV. HOLONOMIC MANIPULATION 

We now proceed to deal with the manipulation of 
the polariton qubit with non-adiabatic holonomic gates. 
Without loss of generality, we set <5 = 0 in our consid¬ 
ered polariton qubit. When a transmon qubit is driven 
by suitable microwave fields, single-qubit quantum gates 
for the defined polariton qubit can be induced. Consid¬ 
ering that a transmon is driven by a microwave field as 
Hd = y/2f(t)a x . It is shown that choosing suitable fre¬ 
quencies in /(f), one can control the transitions between 
between different eigenstates [29] . Here, we are interested 
in the three lowest eigenstates, to keep only the couplings 
between states that form the V-configuration interaction, 
the driven field is chosen as 

fit) = fix cos(wit) + fi 2 cos(w 2 t + ip), (9) 

where f l n is the amplitude of the nth component in the 
driven field with the frequency of ui n , and ip is a pre¬ 
scribed phase difference. Then, we write the full Hamil¬ 
tonian Hi = H,ic + Hd in the subspace spanned by {|G), 
|—), |+)}. Note that this subspace is the eigenspace 
of Hjc, and thus it will be in a diagonalize matrix 
with the elements being its eigenvalues. Meanwhile, as 
(±|(Ta;|G) = (G|cr x |±) = ±l/y/2, the matrix elements of 
Hd will be (±\Hd\G) = (G|ifd|±) = ±/(t). Therefore, 
the full Hamiltonian can be written as 

( E 0 -f(t) fit) \ 

Hi = -fit) £i,_ 0 . (10) 

V fit) 0 E lt+ ) 

In the interaction picture with respect to Hjc, the 
driven fields can be written as 

H[ d = ( -fity^--* o o 

V fit)e^t o 0 
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Setting ui\ = u>o,~ and oj 2 = wothe above Hamiltonian 
reduces to 


( 0 -D 1 D 2 

0 0 

\ D* 0 0 


where D\ ss Ui + VL 2 e l ^ 29t+v ' ) ~ fli and D 2 ~ flie~ 2l9t + 
^l 2 e lip « n 2 e llp . The two step approximations are hold 
under the justification of the rotating wave approxima¬ 
tion. The first one corresponds the omission of the os¬ 
cillating terms with much larger frequencies (~ 2uj\, 2 
2 g). The second one will be hold when g (fii,^). 
Therefore, the driven dynamics of the polariton qubit 
will be governed by 

/ 0 -fii n 2 e iv \ 

H ld = - -fii 0 0 , (11) 

1 \ n 2 e~^ 0 0 / 


as shown in Fig. lb. 

In order to investigate the qubit dynamics under the 
above Hamiltonian, we renormalize it to 

Hlv = | (sin d -e^\G)(+\ - cosf|G)<-| + H.c.) (12) 



4 t/(2w) 



where the effective Rabi frequency is £ = -y/fij-j- f2|, and 
tan(0/2) = tt 2 /S}i. The eigenstates of the Hamiltonian 

m are ig), 

\b) = sin^e“ lv, |+)-cos^|-), 

\ d ) = cos^|+) +sm^e‘ v \-), (13) 


FIG. 2. Fidelity dynamics of the Hadamard gate as a function 
of £t/(27r) with the initial state of the polariton qubit being 
|+), where the red dashed and blue lines being plotted without 
and with decohenrence, respectively. 


holonomic gate. Therefore, the non-adiabatic holonomic 
single-qubit gate can be obtained as 


and thus its dynamics is governed by 

^ = |(|0)<6|+H.c.), (14) 

which means that, in this eigenbasis, the dynamics of the 
three level system can be viewed as a resonate coupling 
between the states |G) and |6) while decouples from the 
’’dark” state | d), as shown in Fig. lc. 

Therefore, the evolution operator U = 
exp i f* Hi v dt'^J realizes the holonomic gate un¬ 
der certain conditions, the reason of which is explained 
as the following. First, in the qubit subspace, the qubit 
states evolve according to \ip±( T )) = U |±). When the 
duration of the Hamiltonian should satisfies f 0 £dt = 2i r, 
the qubit states undergo cyclic evolution. Second, at any 
time, (ip a (t)\Hi\ipb(t)) =0 (a, 6 G {±}). This justifies 
that {IV’-i-(t)}, meets the parallel-transport 

condition in the qubit states subspace, and thus the 
evolution is purely geometric. Note that this condition 
is met due to the structure of the Hamiltonian instead of 
the slow change of parameters in the adiabatic evolution 
case. Under the above two conditions, projected the 
evolution operator onto the qubit space defines the 


U{9,<p) 


cos 8 sin0e 193 \ 
sin 8e l9> — cos 8 J ’ 


(15) 


where 8 and ip can be controlled by the driven field, and 
thus arbitrary single-qubit gate can be obtained by choos¬ 
ing suitable parameters. For example, 8 = 7t/ 4 and <p = 0 
implements the Hadamard gate. 

The performance of this gate can be evaluated by con¬ 
sidering the influence of dissipation using the quantum 
master equation of 

P = #1] + |^(o) + ^£(cr~) + “~U(<T Z ), ( 16 ) 

where p is the density matrix of the considered system, 
C(A) = 2 ApA^ — A^Ap— pA^A is the Lindblad operator, 
k, Ti and T 2 are the decay rate of the cavity, the decay 
and dephasing rates of the qubits, respectively. We con¬ 
sider the Hadamard gate as a typical example and mod¬ 
ulate £ = 5/20 so that £ <C g is well-satisfied to fulfil the 
requirement of g (fii, Cl 2 ) in the derivation of Hamilto¬ 
nian in Eq. m- For icy cs 2-7T x 8 GHz, the cavity decay 
rate is ac ~ 2-7T x 7 kHz [30| and the capacitive transmon¬ 
cavity coupling strength can be g = w r /20 [3lj| . As 
8 = 7t/4 for the Hadamard gate, fli ~ 2.41U2 — 0.924£. 
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Moreover, for a planar transmon qubit couples to a 
transmission-line resonator can, relaxation and coherence 
times of 44 and 20 /js are reported [HI, which leads to 
Ti ~ 27 t x 8 kHz and T 2 ~ 27r x 3.5 kHz. As k , Id 
and r 2 are all on the same order, for simplicity, we treat 
them equally as Tj = T 2 = k = 27rx8 kHz. Sup¬ 
pose the qubit is initially in the state of |+), we evaluate 
this gate by the fidelity defined by Fh = {4’f\p\tpf) with 
| i/jf) = (|+) + | —))/v / 2 being the ideally final state un¬ 
der Hadamard gate. We have obtained a high fidelity of 
Fh — 99.5% at £t/(27r) = 1, as shown in Fig. 2, where 
the red dashed and blue lines are plotted without and 
with decohenrence, respectively. Note that to justify our 
theoretical treatment, we faithfully simulate the quan¬ 
tum process using the original full Hamiltonian Hi in 
Eq. (11011 . and thus do not introduce any approximation. 


V. CONCLUSION 

In summary, we propose a scheme for quantum com¬ 
putation with polariton qubit in circuit QED using non- 
adiabatic holonomy, which is inherently fast and robust. 
In particularly, the polariton qubit is shown to be robust 
against arbitrary low-frequency noise due to its near sym¬ 
metric spectrum, which can also be convenient manipu¬ 
lated by external microwave driven fields in a holonomic 
way. Therefore, our scheme presents a promising geomet¬ 
ric way of manipulating polariton qubits for solid-state 
quantum computation. 
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